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ABSTRACT 

■ We have re-examined the nature of the cluster galaxy luminosity function using the 
| data from the Edinburgh-Durham Southern Galaxy Catalogue and the Edinburgh- 

Milano Redshift Survey. We derive a best fit luminosity function over the range 

■ —18 to —21 in M(bj), for a composite sample of 22 of the richer clusters that has 
O ! M(bj)* = -20.16 ± 0.02 and a = -1.22 ± 0.04. The dominant error in these values 
C****" ■ results from the choice of background subtraction method. From extensive simulations 

we can show that when the LF is fitted over this narrow range, it is difficult to dis- 
r-j . criminate against bright values of M* in the single cluster fits, but that faint values 

O provide a strong test of the universality of the luminosity function. We find that all the 

individual cluster data are well fit by a Schechter function with a fixed at —1.25, and 
that <10% of these have fitted values of M* that disagree from the average at the 99% 
confidence level. We further show that fitting only a single parameter Schechter func- 
tion to composite subsets of the data can give erroneous results for the derived M* , as 
might be expected from the known tight correlation between M* and a. By consider- 
ing two parameter fits, the results of Monte-Carlo simulations and direct two-sample 
X 2 tests we conclude that there is only weak evidence for differences between the data 

■ when broken down into subsets based on physical properties (Bautz-Morgan class, 

richness, velocity dispersion): from our simulations, only the evidence for a difference 
between subsets based on velocity dispersion may in fact be significant. However, we 
find no evidence at all that a Schechter function is not a good model for the intrinsic 
cluster luminosity function over this absolute magnitude range. Models that invoke 
strong evolution of galaxy luminosity of all galaxies within clusters are inconsistent 
with our results. 
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1 INTRODUCTION 

The properties of the cluster galaxy luminosity function 
(LF) are still a matter of active debate. On the one hand 
considerable stress has been laid on a single analytic expres- 
sion for the LF, as determined by Schechter (1976), which it 
is argued holds true for all clusters and the field (e.g. Abell 
1975). On the other hand, much work has been concerned to 
establish differences in the LF as a function of cluster mor- 
phology. Indeed, a universal LF is hard to explain theoreti- 



cally. For example, according to the cannibalism line of ar- 
gument (e.g. Hausman and Ostriker 1978), cD galaxies grow 
at the expense of other fainter galaxies - thus such a model 
predicts a deficit of bright galaxies in the most dynamically 
evolved clusters. The brightness of the 1st ranked cluster 
galaxy is thus related to the dynamical state of the clus- 
ter and the LF is expected to correlate with Bautz-Morgan 
(B-M) type, as this measures the contrast of the 1st ranked 
galaxy compared to the other bright galaxies. Such a deficit 
in bright galaxies is also expected if tidal stripping in the 
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outer halos of galaxies, due to interactions, is taking place 
(Richstone 1975, 1976; Merritt 1983; Malumuth & Richstone 
1984). Only in the special case in which the luminosity func- 
tion is determined during the formation period of the cluster, 
with little further dynamical evolution taking place, are no 
correlations of the LF with other cluster properties expected 
(Merritt 1984, 1985). 

The most popular approach taken to the study of the 
LF is to fit Schechter's (1976) parameterized form, given by 

n(L)dL = n*{L/L*) a exp{-L/L*)d{L/L*), 

where n* is a normalization factor, a is the power-law slope 
at the faint end, and L* (or equivalently M*) is the charac- 
teristic luminosity (magnitude). Schechter(1976) determined 
the best-fit parameter values to be M* — —21.41 (J band) 
and a = —1.25. In subsequent studies there has been claim 
and counter claim regarding the properties of this function. 
In a study of 12 very rich clusters, Dressier (1978) concluded 
that M* does vary significantly from cluster to cluster and 
that it is correlated with the absolute magnitude of the first- 
ranked galaxy, in the sense consistent with the trend ex- 
pected in the cannibalism model. Lugger (1986) studied a 
sample of 9 Abell clusters and concluded that while the mean 
values of M* and a are in good agreement with those of 
Dressler(1978) and Schechter(1976), there was no evidence 
for a correlation of M* with cluster morphology, central den- 
sity or magnitude of the first-ranked galaxy. The examina- 
tion of a further 14 rich clusters by Colless (1989), draws 
very much the same conclusions as Lugger(1986) - with no 
variation in the cluster LFs with B-M type or richness seen. 
Simulations suggest that variations of more than 0.4 mag in 
M* or 0.15 in a are ruled out by these data. Colless did find 
a marginally significant result such that the LF of the high 
velocity dispersion clusters has fainter M* compared with 
the low velocity cases. 

There are potential criticisms of these previous studies: 
first, as pointed out by Lugger (1986), the early studies such 
as that of Dressier (1978) did not use a consistent cluster 
radius or limiting absolute magnitude in comparing differ- 
ent cluster LFs. Secondly, the best previous studies (those of 
Dressier 1978, Lugger 1986 and Colless 1989) confine their 
selection to Abell clusters, and typically only the richer ex- 
amples of these systems. Since there are potential system- 
atic biases in the way in which Abell clusters were selected 
(for example, cD dominated clusters are much easier to de- 
tect by eye), these studies cannot necessarily be held to be a 
true reflection of the properties of all clusters. Last, previous 
studies have relied on global number counts to correct for 
the non-cluster background present, and have used counts 
often taken in different passbands and from different sur- 
veys. This can lead to systematic errors due to differences 
in the calibrations of the data the LF is derived from and 
the data the counts are taken from. 

Here we have re-examined the LF afresh using a sample 
of 46 clusters, selected automatically from digitized data. As 
a result this study has certain advantages over previous work 
in the field. The clusters are selected from photometrically 
calibrated digitized scans of photographic plates using an 
automatic search algorithm. Hence our data should be more 
uniform than the Abell catalogue. The sample covers a wide 
range in richness and many of the clusters have well deter- 



mined velocity dispersions, allowing us to test for variations 
in the LF as a function of intrinsic cluster properties. 



2 THE CLUSTER SAMPLE 

2.1 The catalogue and the redshift survey 

The sample of clusters used in this paper is drawn from the 
Edinburgh/Durham Cluster Catalogue (hereafter EDCC: 
Lumsden et al. 1992), which in its turn is derived from the 
Edinburgh/Durham Southern Galaxy Catalogue (hereafter 
EDSGC: eg Heydon-Dumbleton et al. 1989). The EDSGC 
consists of digitized COSMOS scans of 60 UK Schmidt Tele- 
scope Ill-aJ survey plates centred on the south galactic pole, 
and is nominally 95% complete to bj = 20 with less than 
10% contamination. A cluster sample was derived from the 
galaxy data in the EDSGC as outlined in Lumsden et al. 
There are over 700 potential clusters in the EDCC, and this 
sample is the first published machine based cluster sample. 
The objective nature of the construction of the EDCC im- 
plies that it does not suffer from the systematic subjective 
biases that appear in the southern extension of the Abell 
catalogue covering the same area (Abell, Corwin & Olowin, 
1989). The EDCC is complete within its selection criteria to 
a limiting magnitude for the tenth ranked cluster member 
of bj = 18.75. 

For a smaller sub-sample of the EDCC clusters, we ob- 
tained redshifts with the multi-slit spectrograph EFOSC on 
the ESO 3.6m and with the multi-fibre spectrograph Aut- 
ofib at the 3.9m Anglo-Australian Telescope. The nature of 
this sample is discussed in Nichol et al. (1992) and Guzzo 
et al. (1992) and the redshifts are published in Collins et 
al. (1995). Approximately 10% of those clusters observed 
proved to have significant contamination from other groups 
along the line of sight at significantly different redshifts. We 
have excluded those clusters with 'projection effects' as de- 
fined in Collins et al. (1995) from our sample. We note that 
the remaining clusters do not have a well defined cut-off in 
cluster richness (unlike the sample of Nichol et al), since we 
have included all the clusters for which we actually obtained 
redshifts. 

From this dataset, we have excluded those clusters 
where the cluster background radius (as defined in section 
3.3) overlaps the boundaries of the EDSGC for ease of com- 
putation, and those clusters where their respective Abell 
radii are entirely superposed, since there is then no well de- 
fined cluster magnitude distribution. We also individually 
checked all the other clusters where the two clusters' Abell 
radii overlapped since the deblending procedure described 
in Lumsden et al. can occasionally cause the derived cluster 
centroid to move away from the 'true' position: we removed 
those clusters where the catalogue positions differed from 
both the centre derived from a visual inspection of the plate 
and the centroid of those galaxies for which we have red- 
shifts by a significant fraction (~30%) of the Abell radius. 
For blended clusters that pass this test, we assign galaxy 
membership to the individual clusters using the method de- 
scribed in Lumsden et al. (1992). Lastly, to be confident 
that we have removed all clusters that may have contami- 
nating projection effects, we choose to use only those with 
a minimum of 4 cluster redshifts, and where the number of 
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galaxies with redshifts outside of the cluster is less than the 
actual number of bona-fide cluster members. 

The parameters of the remaining 46 clusters (as taken 
from Lumsden et al. 1992 and Collins et al. 1995) are given 
in Table 1. For all of the clusters, we have considered cluster 
members to be those galaxies within one Abell radius (where 
we have taken ta = 1.5/i _1 Mpc and we use the value for 
the Hubble constant h — 1 throughout *). Our sample is 
the largest single database of clusters which has been used 
to study the LF to date. 



plate corners) will consequently affect the amount of satu- 
ration present. We do not have sufficient CCD sequences to 
map these effects fully however, although we note that the 
overall photographic intensity-magnitude relation has been 
corrected for these effects (see Heydon-Dumbleton 1989). 
These effects may in part give rise to the observed scatter in 
Fig. 1. We expect the magnitude of any error caused by us- 
ing this average calibration to be small, since the deviations 
from the fit are small. 



2.2 Magnitude calibration 

The EDSGC was calibrated from photometric CCD se- 
quences assuming that the plate and CCD magnitudes were 
linearly related (Heydon-Dumbleton et al. 1989). However, 
for bright objects whose profiles are either strongly peaked 
(such as elliptical galaxies) or stellar, there is a regime in 
which the measured COSMOS magnitude is less than the 
true magnitude, caused by lack of dynamic range within the 
measuring machine itself. This 'saturation' occurs for ellip- 
tical galaxies at bj ~ 17 (note that the plates themselves 
only saturate at much brighter magnitudes). Fainter than 
this magnitude, and for most late type galaxies, the original 
calibration is correct, so that previously published results 
from the EDSGC are unchanged by these corrections. How- 
ever, since most clusters are dominated by elliptical galaxies 
at the bright end of the LF it is vital that we correct for this 
effect or else we will induce errors in the measured value of 
the 'break' in the LF for the nearby clusters. 

We have calculated a correction for this effect, by com- 
paring the catalogued EDSGC magnitude against the ob- 
served magnitude for 451 galaxies present in the CCD se- 
quences used to calibrate the EDSGC initially (Heydon- 
Dumbleton et al). Since for faint magnitudes the slope 
should be linear, but at bright magnitudes it should diverge 
from this, we used a quadratic fit to the data. From this we 
derive a fit with the form (see Fig. 1) 

bj (EDSGC) — bj(CCD) = (1.06 ±0.26) x 10"% (EDSGC) 2 
-(0.46 ± 0.09)63 (EDSGC) + (4-96 ± 0.90). 

This fit should be reliable in the range 15 < bj < 21. For 
most of this range, the linear component present is actu- 
ally the most important. This slope has been commented on 
previously by Heydon-Dumbleton (1989), and gives rise to 
the observed slope in the fit at fainter magnitudes, and is a 
result of the definition of measured photographic magnitude 
used in the EDSGC. As can be seen the effect is small (the 
overall effect on the composite LF described in section 5.1 
for example is 0.2 magnitudes). However, for the purposes 
of the LF, we include all terms, so that our magnitude scale 
is more closely matched to the CCD based magnitude scale. 

It should be noted that this correction is an average, 
and may not represent every cluster absolutely. Since the 
saturation effect is a function of the photographic density, 
anything else that modifies this (such as vignetting, emul- 
sion flaws, plate fogging or emulsion desensitisation near the 

* Where h =Ho/100, and Ho is the Hubble constant in units of 
km _1 s _1 Mpc _1 



3 THE CLUSTER GALAXY LUMINOSITY 
FUNCTION 

3.1 Form of the luminosity function 

As is common with studies of the cluster galaxy luminosity 
function we use the standard Schechter function (Schechter 
1976) as the analytical model to our data. This function 
has the following form when expressed in terms of absolute 
magnitude: 

n(M)dM = knx 

exp{fc(a + 1)(M* - M) - exp(fc(M* - M))}dM, (1) 

where k = lnlO/2.5. Here M is the absolute magnitude, 
derived from the measured apparent magnitude in using, 

M — m — 42.384 + 5 log h - 5 log z - K(z) - A(bj), 

where K(z) — 4.14,2 — 0.44z 2 is the K-correction suitable for 
the bj passband for the mix of galaxies typically observed in 
low redshift clusters (Ellis 1983), and A(bj) is the effective 
extinction in this passband. This latter quantity is derived 
using techniques described in Nichol & Collins (1993) for 
each cluster (see also appendix A). In common with other 
studies of the luminosity function we assume that go = 1 in 
deriving the distance modulus used, and we have adopted 
h = 1 for simplicity as well. The LF normalisation n* is con- 
strained using the actual observed number of cluster galax- 
ies, and we carry out either a two parameter fit to M* (the 
absolute magnitude of the 'break' in the LF) and a (the ex- 
ponent of the power-law slope at fainter magnitudes) , or a 
one parameter fit to M* with a fixed where required. 

3.2 The Composite Luminosity Function 

There are too few galaxies per bin in the individual LFs to 
fit both M* and a simultaneously. Since we wish to search 
for deviations from some 'mean' LF, it is valuable to con- 
sider coadding all the LF data into an average or composite 
luminosity function, in a way that allows us to estimate both 
parameters simultaneously. The easiest manner in which to 
do this is to consider the composite luminosity function de- 
fined by Colless (1989) to be 



where Ri is the 'richness' of the ith cluster, n c ius,j is the 
number of clusters that actually contribute to the jih bin 
of the composite LF, Nij is the background corrected num- 
ber of galaxies in the jth bin of the zth cluster's LF and 
R c — Yi c i us Ri. We use a different definition of Ri to that of 
Colless, since he used the total number of galaxies brighter 
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than M 



-19 and we use the background corrected 



number of cluster galaxies between M{bj) = —19.5 and 
M(bj) = —20.5. This is because our sample contains more 
distant clusters on average than the Colless sample, and we 
would otherwise discard a significant fraction of the data as 
being objects where M = — 19 was beyond our chosen com- 
pleteness limit. We also impose an upper limit in absolute 
magnitude when calculating the richness to avoid the prob- 
lems with very bright galaxies discussed above. For typical 
values for the LF (M* = —20 and a — —1.25), the relation- 
ship between our definition of richness and that of Colless is 
Jii(CoUess)~ 2R t (here). 

Similarly, the errors for the composite LF are calculated 
according to 



SN cj = 



Rc 



1/2 



(3) 



where the SNij are derived according to the error given in 
Section 3.4. 



3.3 Background galaxies 

In order to derive the LF, we must be able to define the 
expected contribution from background galaxies, nt ac k, ac- 
curately. We have considered three approaches to this prob- 
lem. First, we used the counts in an annulus, centred on the 
cluster, with an inner radius of one degree and a width of 
one degree. This fixed radius was adopted for simplicity - 
using a fixed metric radius as opposed to a fixed angular 
radius gives essentially the same results. This method has 
the advantage that the correction is strictly local, and hence 
the variation in field counts seen in different regions of the 
EDSGC is correctly accounted for. For the other two meth- 
ods we used the the global counts as derived from the whole 
of the EDSGC, and the average of the counts from the annuli 
from the first method described above. The first of these is 
similar to the procedure used in previous studies, but since 
our global counts actually come from the same data the LFs 
themselves are constructed from, we remove the possibility 
of errors arising from discrepancies in the calibration of the 
LF and the global counts. The second has the potential ad- 
vantage that it excludes the rich clusters from the counts 
and hence is a better measure of the 'true' background. In 
Fig. 2 we show the derived counts from the catalogue as a 
whole together with the average of the counts derived from 
the annuli around the clusters. There is essentially little dif- 
ference between these on average. The difference that does 
exist is at the level expected from the presence of real large- 
scale structure. There is little sign that excluding the rich 
clusters cores makes the background as calculated from the 
annuli lower than the global counts. 

Previous estimates of the background count for use in 
the LF have been derived from either separate 'eye-ball' 
galaxy catalogues (e.g. Dressier 1978, Lugger 1986) or have 
been based on data from only very small areas (Colless 
1989). We note in passing that the slope of the number 
counts that we find is larger than that assumed in previ- 
ous studies of the LF. This is however consistent with Mad- 
dox et al. (1990), the only other measure of the counts over 
significantly large areas of the sky that is derived using a 



measuring machine. The fit given by Colless (1989) appears 
to overcount the true background by ~ 20%. 



Lastly, we note that the derived number counts shown 
in Fig. 2 drop off more rapidly than m 0,45 at bj > 20. This 
disagrees with other observations (eg the deep CCD counts 
presented in Metcalfe et al. 1991), and indicates that the 
overall completeness limit of the EDSGC is bj ~ 20. How- 
ever, inspection of the number counts on individual Schmidt 
plates also shows that this limit varies from plate to plate. 
We have allowed for a possible variable completeness with 
plate in our work by examining the number counts on each 
plate separately and assigning completeness limits relative 
to the 7 deepest fields in the survey. From this we find that 
the best plates are complete to bj = 21 and the worst to 
bj = 19.5. We do not however use data fainter than bj = 20 
from any field. For those fields where the adopted complete- 
ness limit is brighter than this, we do not use data fainter 
than that limit. In calculating the background correction 
required when using the global number counts we have in- 
terpolated the actual measured data, rather than fitting a 
single power law. The latter is a poor representation to the 
actual data, since the true number counts 'turn-over' near 
bj ~ 19 (eg Metcalfe et al.). This regime is important for our 
work since for clusters with z ~ 0.1, bj • = 20 is equivalent to 
an absolute magnitude of ~ —18. 



We also have to derive suitable error estimates for these 
background counts. There are two sources of error, which are 
present in all the methods given above. The first is simply 
the Poisson uncertainty in the counts and can be represented 
simply as ^friback- The second is systematic and represents 
the presence of large-scale clustering within the sample. For 
the global correction method, although we would assume 
that the slope of the background number counts is constant 
over the survey region, the normalisation of these counts 
changes with position. Colless (1989) and Lugger (1986) fol- 
low Geller & Beers (1982) in assuming that the background 
is only known to 50%. For the background derived from the 
local annuli we can estimate the error rather more accu- 
rately however. If we take all of these annuli, and find the 
scatter about the mean of the number counts, then we can 
derive the systematic error for this method. Since all the an- 
nuli have the same area any number density deviations are 
purely due to the presence of large-scale structure. The de- 
rived variance is also shown on Figure 2. Comparison of the 
expected contribution from Poisson error and the total er- 
ror observed indicates that the mean field-to-field variance is 
25% for bj > 17. The largest excursions seen are <50% over 
the same magnitude range. Therefore, using the 50% vari- 
ation in what follows, although conservative, does account 
for all the observed variations, whereas using the 25% mean 
variation is a better representation of the data as a whole 
but will underestimate the error in the background for a few 
clusters. We demonstrate in section 5.1 the effect of using 
this smaller error on the field counts compared to the value 
favoured by Geller & Beers. We note that Dressier (1978) 
also found that a mean variation of 25% was a good fit to his 
data, and this result is compatible with the known angular 
clustering in the EDSGC (Collins, Nichol & Lumsden 1992). 
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3.4 Fitting Methods 

We have considered only a x 2 squared fitting method in 
deriving the parameters M* and a. If we let n c ; us be the 
difference between the total counts in a bin and the expected 
background in the bin (however the latter is calculated), so 
that 

then we seek to minimise the function 
2 _ ^ (ridus - rifit) 2 . . 

where, 

n"(M)AM 3 
n/« = n(M)AM + ^ 

and 

cr(nfit) = [o-(n t0 i) 2 + cr(rti, acfe ) 2 ] 1/2 . 

Throughout we use a bin width, AM, of 0.5 magnitudes. 
The presence of the second derivative in the form of Ufa 
takes account of the finite bin width used, since the number 
counts change extremely rapidly (Schechter 1976). The first 
term in the error represents the Poissonian error in the total 
observed counts. Therefore, a(ntot) = n]^ 2 . The second term 
must take account of both Poissonian error and systematic 
errors in the background as noted in section 3.3. This has the 
form a(n back ) = max(n^ 2 fe , xn back ), where x = 0.5 (except 
where explicitly noted that x = 0.25). 

Although for many of the clusters data exists with M 
brighter than —21 and fainter than —18, there are good rea- 
sons for not using such data in deriving values of M* and 
a. First, for the fainter galaxies, the correction for the back- 
ground becomes more difficult since the background fraction 
per bin rises steadily. For any given cluster we impose two 
limits on this: that any data with m(bj) > 20 be discarded, 
for reasons outlined in section 3.3, and that any bin for 
which the number of background galaxies exceeds the clus- 
ter galaxies also be discarded. In practice, the former limit 
essentially ensures the latter for faint galaxies in any event. 
For bright galaxies there are also problems we wish to avoid. 
It is desirable, as has been noted previously by, e.g. Lugger, 
to exclude the brightest cluster galaxy when deriving the 
LF. We automatically do this since we discard all galaxies 
brighter than M — — 21 (typically of order five galaxies in 
the richer clusters). There are more specific reasons for do- 
ing this in our case as well. First, the saturation correction 
is poorly defined brighter than bj = 15, and hence the mag- 
nitudes of the galaxies brighter than this are likely to have 
larger errors. At z — 0.05, this magnitude corresponds to 
M = -21.5, and at z = 0.1 M = -23. Therefore, it is clear 
that galaxies brighter than M = —21 should not be used be- 
cause their photometry is suspect. Secondly, COSMOS data 
is deblended where object mergers occur (cf the discussion 
in Lumsden et al. 1992), but the accuracy of the magnitudes 
of the deblended data is less reliable than for single isolated 
objects. The brightest galaxies in clusters are almost always 
blended with fainter satellite systems in the EDSGC data, 
so again the photometry of these objects is not reliable. Dis- 
carding the brightest bins of the LF guards against problems 
of this kind. It is possible that similar problems may also be 
present in previous derivations of the LF, since the work 



described in the introduction also relied extensively on pho- 
tographic plate material. In any event we recommend that 
results which rely on accurate magnitudes for the brightest 
galaxies in any LF derived from photographic plates should 
be treated with caution. We have also placed constraints on 
the values of M* and a that the fitting routine can accept, to 
prevent any possibility of the fitting process 'running away.' 
These limits are that M* should lie between —18 and —26 
and that a should lie between and —2.5. Therefore quoted 
values for the fits that are near these limits are indicative 
of a very shallow minimum in the fit, and hence that the 
reliability of the fit is likely to be low (whatever the quoted 
formal probability that the fit is good) . Lastly, we note that 
we sum bins to ensure that there are > 5 galaxies per bin 
in the data being fitted, as is standard for x 2 minimisation 
techniques. 



4 RELIABILITY OF THE FITTING 
PROCEDURE 

Before we can consider the actual data, we need some esti- 
mate of the reliability of our fitting procedure and the true 
error on these fits. To achieve this, we conducted extensive 
Monte-Carlo simulations. In particular, we carried out three 
sets of simulations: first, to estimate the effect that errors 
in the assumed background had on the fit parameters as a 
function of redshift; secondly, to derive a measure of the true 
error distribution for a one parameter fit to a model LF as a 
function of richness, for Ri typical of the single cluster values 
we found; lastly, to derive the error distribution for the two 
parameter fits to the composite sets. These latter two will 
then allow us to estimate the true likelihood that observed 
departures from the composite LF are in fact truly signif- 
icant. In carrying out these simulations, we have assumed 
that M* = -20.2 and a = -1.25 for the model of the LF 
(in accordance with the results derived in section 5.1). 

In order to check the effect of the background subtrac- 
tion alone, we initially assumed that the form of the LF 
was known perfectly, and only the background was allowed 
to vary. We used the global background as derived above, 
and rescaled it at random so that the variance was 25% as 
was found to be the case in practise in section 3.3. Thus 
after normal background correction, the random LF would 
have either too many or too few background galaxies sub- 
tracted. By deriving fit parameters for many such datasets 
we can derive the average error induced in those parameters 
as a function of cluster richness and redshift. In general, for 
z < 0.05, the effect of errors in the background on the fit 
parameters is negligible, even for poor clusters. At higher 
redshifts, M* can be measured to an accuracy of 0.05 mag- 
nitudes if Ri > 30 (ie the average of the measured M* from 
the simulations is different by this amount from the assumed 
input M*). At z > 0.05 and Ri < 20, the scatter in the de- 
rived value of M* becomes larger (~ 0.5 — 1.0 magnitudes). 
The values derived here should be the dominant source of 
error when the richness is very large (since then the LF 
is known almost perfectly). Since the composite luminosity 
functions given all have effective values of the richness that 
are much greater than 100, it is clear that the empirically 
derived errors for the effect of different background correc- 
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tions given in section 5.1 are in good agreement with the 
results of these Monte-Carlo simulations. 

More generally, we need to map out true estimates of 
the likely scatter in the fits to M* when the LF was drawn 
at random from the same model LF used above. In this case, 
since we draw individual galaxies at random from the par- 
ent LF, we also include random fluctuations in the LF itself. 
For low to moderate richness clusters, it is these fluctuations 
that dominate the scatter in the fits to M* rather than the 
background correction method. Although the resultant val- 
ues are therefore independent of redshift, we have adopted 
2 = 0.1 as typical of our clusters, in order to include a suit- 
able background contribution. This background is derived 
in the same manner as described above. Table 2 gives the 
limits outside of which 90%, 31.5%, 10%, 5% and 1% of the 
data lie. One of the most noticeable features of these sim- 
ulations is that the distribution of fitted values for M* is 
skew, with a strong tail towards brighter magnitudes. An 
example of this is shown in Figure 3, where we present the 
distribution of fits to M* for two richnesses, Ri = 10 and 
Ri — 50. This trend is due to the restricted range of abso- 
lute magnitudes that are actually used to fit against. Since 
we set M* = —20.2 in accordance with the result we find 
in section 5.1, there are only two bins in the fitted LF that 
actually sample the break. At low richness values there are 
often insufficient counts per bin to truly define this break, 
and hence, since there is no constraint from the actual data, 
the fitted value of M* can take quite large negative values 
without changing the \ 2 of the fit significantly. In practice, 
the limit for M* at the bright end is set by our constraint 
that M* > —26. Clearly, for the richer sample shown in 
Figure 3 the trend towards a strongly skewed distribution is 
much less evident as expected. This skewness also shows in 
the mean M* derived from the simulations, with this value 
being the input value of —20.2 at large richness but diverging 
to —20.8 at Ri = 10. This same behaviour is also demon- 
strated (in a slightly different fashion) by the simulations of 
Colless (1989), where he shows that the probability that a 
two-sample \ 2 lest can reject the possibility that two clus- 
ters are drawn from the same underlying population is also 
skew with respect to M* . We can further confirm the reason 
for this skewed distribution, by looking at the distribution 
of M* when the absolute magnitude range being fitted is 
—23 to —18. The results for Ri = 25 are given in Table 2 
and also shown in Figure 3. It is easy to see the effect that 
including the brighter galaxies has on fixing the location of 
the break in the LF. The slight skewness left in the distribu- 
tion now reflects the fact that the total number of galaxies 
with M < —21 in these low richness systems is rather few. 
The derived error bounds including these bright galaxies are 
most similar to those for a much higher value of Ri (~ 40) 
when fitted over the normal range in M. The key result 
of these simulations therefore is that it is difficult to rule 
against overly bright values of M* as being from the same 
distribution as the composite. However, those clusters for 
which M* is rather fainter than the composite are strongly 
discriminated against. We note that the results of these sim- 
ulations can also be applied to the data of Colless (1989), 
and that the simulations of Dressier (1978) are not applica- 
ble to the case of fitting over a restricted magnitude range 
(or rather that the Dressier results would be closer to those 
achieved in practice but for a much larger true richness) . 



It is worth also considering how sensitive these results 
are to our assumed value of a. We tested this using a = 1.2 
for Ri ~ 20, typical of our sample. We found that the 
changes are small (< 0.1 in the faint limit on M* , and < 0.2 
in the faint bright limit) and always act to push these lim- 
its brighter. Since the cluster LFs we can exclude are those 
which have faint M* , it is clear that adopting a lower value 
of a will only make this difference larger. The formal aver- 
age of our composites using all methods of calculating them 
give a = —1.22 (Table 3). This implies that the probability 
that any LF agrees with the derived composite as given by 
the simulations is slightly conservative. 

Since we have carried out two parameter fits in section 
5.1 and 5.3, we also need to know the likely joint distribution 
of a and M* for the case of rich systems alone (ie compos- 
ite LFs). We therefore carried out simulations as described 
above but in this case applied a two parameter fit to the LF. 
In this instance the best way to present the data is graphi- 
cally, and the results of the simulations are shown in Figure 
4. The contour levels plotted correspond to the same prob- 
abilities given for the one parameter case in Table 2, except 
that we exclude the contour that contains 10% of the data 
in this instance. The simulations have been binned into 0.05 
in a and 0.1 in M* for presentation. 

It is clear from Figure 4 that M* and a are highly corre- 
lated. In what follows we will present formal error estimates 
on each fit, derived by determining the value(s) at which 
X 2 = Xmin + 1. As will be seen, there is a considerable 
difference between these formal errors and both the error 
distributions given in Table 2 for the one parameter fits and 
the distribution shown in Figure 4 for the two parameter fits. 
This difference largely stems from the correlation between 
M* and a, since the errors in these are also therefore clearly 
correlated. These formal errors are included for complete- 
ness only (since they do provide a measure of 'goodness- 
of-fit'). We will only use the Monte-Carlo simulations and 
two-sample x 2 tests when assessing the true deviations that 
may be present between any given LF and a universal com- 
posite. 



5 RESULTS 

5.1 The Composite Luminosity Function 

The first LF we consider is the composite of all the clusters 
which have Ri > 20. There are 23 clusters in this sample, 
and the richness of the composite was found to be R c ~ 800 
(n*n c i us ~ 1600). This allows us to compare our results 
with previous samples which generally dealt with rich Abell 
clusters. 

For this composite we have considered the effect of us- 
ing the different background correction methods outlined in 
Section 3.3, as well as for the different assumed errors in 
that background subtraction. The results of the fitting pro- 
cess for these different background corrections are given in 
Table 3 and shown in Figure 5. 

The first three rows were derived by fitting both M* 
and a to this sample over the magnitude range (—21 < m < 
— 18) with an assumed field-to-field error in the background 
counts of 50%. The first entry gives the results for the lo- 
cal correction method, the second for the global correction 
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using all of the counts and the third the global correction 
but using only the average of the annuli around the clus- 
ters. Clearly the two global background corrections we have 
used give essentially the same result. Directly comparing the 
data derived from the local and global background correc- 
tion methods using a 2-sample \ 2 test also shows that these 
composites are in good agreement for the absolute magni- 
tude range given above (^99% probability they are same). 
Moreover, since we have not renormalised the LFs before 
comparing using this test, this also indicates that differences 
that do exist between the background subtraction methods 
for individual clusters (as will be shown in Section 5.2) tend 
to average out when combined in the composite. 

The last three rows in Table 3 give the equivalent results 
when the field-to-field background error is assumed to be 
25% as noted in Section 3.3. This change therefore gives the 
same results as before within the derived errors on M* and 
a for the global correction method, and only very slightly 
different in the formal sense for the local correction method. 
This shows that the dominant error in deriving estimates 
of M* or a is actually the Poissonian contribution (cf the 
discussion after equation 4), and that the variation due to 
assuming a different error on the background term alone is 
small. Indeed, the magnitude of these changes found here is 
negligible compared to the expected measurement errors as 
shown in section 4. Henceforth, we consider only the 50% 
field-to-field variation, since it adequately represents all of 
the clusters in the present sample. 

Lastly, it is worth noting the effect that changing 
the magnitude range over which the fits are made has on 
these conclusions. Including the brighter galaxies pushes M* 
brighter. Including the fainter galaxies makes a slight differ- 
ence to the slope but little else. It should be noted how- 
ever that only one cluster actually contributes for M fainter 
than —18, whereas all 23 clusters contribute at M brighter 
than —21. The result derived when the fit is over the range 
—22 < M < —17.5, using the local background correction, 
is also given in Table 3. The low probability that the fit is 
good is due to the fact that a Schechter function is now a 
poor fit to the data at both bright and faint magnitudes. 
We note that a similar (though larger) increase in brightness 
of M* was found by Lugger when she considered mean lumi- 
nosity functions with the brightest cluster galaxies included. 
However, we repeat our caution of section 2.2 that data on 
such bright galaxies (m < 15) cannot be relied upon abso- 
lutely because of problems of saturation and image blending 
(cf the discussion in Colless 1989). It is also not clear there- 
fore whether or not the greater variance from the mean seen 
in the brightest bins, as compared to those in the more re- 
stricted absolute magnitude range over which we normally 
determine the LF, is actually real or merely an artifact of 
the COSMOS measurements. 

5.2 Individual Cluster Luminosity Functions 

We also fitted each cluster LF separately to the Schechter 
function, but this time holding a constant (since there are in- 
sufficient points in many of the clusters to enable more than 
a one parameter fit to be made successfully). We adopted 
a = —1.25 as an appropriate compromise, based both on our 
own results and those previously published. Table 4 shows 
the results of these tests for all 46 clusters that were found to 



have positive values of Ri and Figure 6 shows the actual data 
for those with Ri > 15. We quote results derived using both 
the global and local background correction methods. We also 
tested for differences arising from the choice of background 
subtraction method by comparing the LF distributions di- 
rectly. We used a two-sample \ 2 test to check for differences 
between the two sets of LFs, without rescaling the LFs for 
the possibility that the two background subtraction meth- 
ods may have resulted in a different overall normalisation. 
We have allowed for the fact that the adopted magnitude 
completion limit may be brighter than M = —18, so that 
we compare the LFs in the range —21 to the brighter of —18 
or the plate completion limit. These limits are given in Table 
4. The results of the two-sample \ 2 test are given in column 
10 of Table 4. As can be seen the agreement is good, and 
even for the worst case (E748), we cannot rule out the pos- 
sibility that the two samples represent the same underlying 
distribution. 

We also compared the single cluster LFs with a suitably 
normalised model of the composite (assuming the average 
values of M* and a found above), using a simple \ 2 test 
(given in Table 4), and with the results of the Monte-Carlo 
simulations from Table 2. It can be seen that there are sev- 
eral clusters where the LF derived using global background 
subtraction is a poor match to the model, but the LF de- 
rived using the local background is a much better fit. We 
take this as evidence for the fact that assuming an overall 
global background will not be a good match to every cluster. 
Generally the two tests give similar results, though we note 
that the two-sample \ 2 test is considerably more conserva- 
tive at ruling out individual LFs, as expected. However, all 
those LFs for which the two-sample x test gives a proba- 
bility below 20%, are also found to be a poor match to the 
composite from the simulations. There are two exceptions 
to the corollary however: both E438 and E462 have unex- 
ceptional probabilities from the two-sample \ 2 test, but are 
excluded at the 5% level from the simulations. 

If we include only those clusters where the derived M* 
lies outside the 1% confidence level from the simulations (as 
taken from Table 2) for one or other of the background sub- 
traction methods, then only clusters E400, E429 and E524 
show strong evidence for an LF that differs from the compos- 
ite. All three of these clusters also disagree at the 5% level 
when considering both background subtraction methods. At 
the 5% level, there are several additional clusters that dis- 
agree. Leaving aside those where only one of the background 
subtraction methods gives divergent results (E499, E726), 
we are left with E438, E462 and E742. All of these, except 
E400 and E429, have a value of Ri < 15, so we may rea- 
sonably exclude these as being different on the grounds that 
accurate background subtraction is crucial for these clus- 
ters, and, as found in Section 4, errors of 0.5 magnitudes 
in the fitted value of M* are easily accounted for by an in- 
appropriate background correction. This leaves only E400 
and E429. It can be seen from Table 4 that even the direct 
X 2 comparison of their LFs with the composite might lead 
to them being thought discrepant. The LFs of these clus- 
ters are apparently normal (see Figure 6). However, before 
concluding that these clusters do show a significant differ- 
ence, we first consider the possibility that there may be some 
problem with the data. Examination of the data in Collins 
et al. (1995) shows that all of these clusters have good re- 
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liable redshifts, with the positions of the galaxies for which 
redshifts were obtained agreeing well with the cluster cen- 
tres given in Table 1. Further, we note that for the cases 
of E400 and E429, other clusters appear on the same pho- 
tographic plate and do not have discrepant M* , hence the 
photometry of these systems should be good. Therefore, we 
find that there are at least two clusters that do not have the 
same LF as the derived composite. This is in keeping with 
earlier studies by Dressier (1978), Lugger (1986) and Colless 
(1989) who found that ~10% of their clusters similarly did 
not agree with the mean Schechter LF they derived. 

5.3 Dependence of the LF on Cluster Properties 

The major advantage of our sample over previous studies 
is that it is large enough that we can test the dependence 
of the LF on intrinsic cluster properties. In particular, we 
will consider the possible variations with cluster richness (as 
defined by Ri), Bautz- Morgan class, which we have taken 
from the Abell, Corwin & Olowin catalogue, and velocity 
dispersion (taken from Collins et al. 1995). 

For each case we split our total sample into two, and 
used the local background correction method only. When 
comparing different richness samples we used the cuts 10 < 
Ri < 20 and Ri > 20 (which is of course just the composite 
LF derived previously). For Bautz-Morgan class, we consid- 
ered only those of type I and III, and counted intermediate 
types (e.g. type I-II) as type II which we did not use. In this 
way we sought to maximise the differences between the clus- 
ters to test for any effect on the LF. We split the sample into 
two in velocity dispersion using a cut at 700kms _1 . Unfor- 
tunately, most of the very high velocity dispersion clusters 
(those with v > lOOOkms -1 ) are also the most distant, mak- 
ing it difficult to constrain a, and leaving few bins in the 
data to compare the 2 sets using a x test. We therefore 
were forced to use a lower velocity cut-off for our high ve- 
locity dispersion sample than Colless (1989) did, so do not 
sample the same part of parameter space as he did. 

Since we do not have well defined velocity dispersions 
for all clusters, and those which are not in the Abell, Cor- 
win & Olowin catalogue do not have Bautz-Morgan classes, 
the two subsamples do not comprise all of the observed clus- 
ters. The number of clusters comprising each subsample is 
given in Table 5. In addition, for the velocity dispersion and 
Bautz-Morgan samples we have included two richness cuts, 
at Ri > 20 and Ri > 15. This provides some indication of 
the stability of the result on the clusters making up the ac- 
tual subsamples. We give both one and two parameter fits 
to the Schechter LF in Table 5. For the one parameter fits 
we again used a = —1.25. The results from this analysis are 
also shown in Figure 7. We have also compared the LFs for 
these subsamples against the composite derived in Section 
5.1, and against each other, using a two-sample \ 2 test, as 
well as comparing the subsamples directly with the results 
of our simulations (Figure 4). The results of the two-sample 
X 2 test are also given in Table 5. Since n* can vary between 
these different sets of clusters, we rescaled the binned data 
before the comparison, so that the derived R c values from 
the composites agreed. 

There are noticeable differences between the richness 
cuts for both the Bautz-Morgan class III clusters and the 
high velocity dispersion clusters, with the richer subsamples 



having a notably poorer fit to the standard Schechter LF, 
and being discrepant compared to the composite. However, 
only one cluster contributes to the faintest bin in the richer 
subsample, and only one other is added to that when the 
poorer limit is considered. If we exclude this bin, then very 
similar fits are found for both richness limits and the dis- 
crepancy is removed. It is clear from this that variations in 
the LF due to the last bin alone are not significant. We do 
not consider these two subsamples further. 

Once these discrepant sets are removed, it is clear that 
neither Bautz-Morgan class shows any evidence for devia- 
tion from a universal LF, nor do they differ from each other, 
according to any of the tests applied. The same is true for 
the split according to richness alone, where the fits to both 
samples are very similar within the derived errors. The ac- 
tual 'goodness' of the fit for the lowest richness sample is 
poor, but this most likely reflects difficulties in background 
correction at these low richness values as noted previously. 

For the Ri > 15 velocity dispersion samples, there is 
weak evidence that the two samples differ from each other 
(significance ~ 54%). When compared with the composite, 
neither sample shows any evidence for difference when using 
the two-sample \ 2 test. However, when compared with the 
Monte-Carlo simulations, we find that the higher velocity 
dispersion clusters are clearly unlikely to have been drawn 
from the same basic distribution as the composite (at much 
better than the 1% level). We therefore find that there may 
be evidence for some dependence of the LF on velocity dis- 
persion. 

It is worth considering the effect of including the 
brighter galaxies as well in this analysis. As is readily ev- 
ident from Figure 7, there is little difference between the 
rich and poor samples across the entire range of M, and the 
differences for the high and low velocity dispersion samples 
is largely restricted to the fainter galaxies. Considering the 
Ri > 15 Bautz-Morgan samples however, there are clearly 
fewer bright galaxies in the Bautz-Morgan class I sample 
than in the Bautz-Morgan class III sample. This is not a 
normalisation effect since the two samples have essentially 
the same LFs at fainter magnitudes. Again however, we cau- 
tion against overinterpretation of this difference. 

Lastly, we note that all of the velocity dispersion sub- 
sets, and the low richness subset can formally be excluded 
as being drawn from a model of the composite LF with a 
fixed at —1.25 on the basis of the one parameter fits. Small 
deviations in both M* and a mean that the one parameter 
fits are distinct from the composite derived in section 5.1. It 
is clearly dangerous to compare such subsamples (with their 
relatively small intrinsic errors) using such methods. 



6 DISCUSSION 

6.1 Comparison with Previous Results 

First, we compare our derived estimates for the composite 
LF, for all clusters with Ri > 20, with those derived previ- 
ously. Colless (1989) gives the best fits from previous data 
in his Table 2, scaled into the colour and value of Ho that 
we are using. Excluding Colless' own data this shows that 
the quoted values of M* and a are —19.9 and —1.24 respec- 
tively, with typical errors in M* of 0.5, when the brightest 
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cluster galaxies are excluded from the fit. Colless himself 
finds M* = —20.04 and a = —1.21 for a two parameter fit 
similar to ours. These values are in good agreement within 
the errors with those derived by us. It is clear that the lumi- 
nosity functions of moderately rich clusters, when considered 
across this relatively narrow range in absolute magnitude, 
are indeed very similar on average. 

In addition we can compare our one parameter fits for 3 
clusters with those given by Colless (1989). EDCC clusters 
394, 400 and 124 correspond to his clusters C02, C03 and 
C52. From Table 4 we can see that when a comparison is 
made of similar fitting procedures (ie, x 2 fitting and using 
a global background correction), the agreement is excellent. 
Colless quotes values for M* of -19.76±0.21, -19.86±0.23 
and -20.05 ± 0.33 for C02, C03 and C52 respectively. Al- 
though the original photographic plates used in both studies 
are the same, a different measuring machine was used (the 
APM), with consequent differences in the reduction from 
measuring machine magnitudes to final derived astronom- 
ical magnitudes. Therefore, this agreement for the specific 
clusters that are common to the two studies is highly en- 
couraging. It is worth making some specific comments on 
these clusters however. E394 was found by Colless to be 
discrepant with both the fit by a Schechter LF and with 
a comparison with his composite LF. We note that Table 
4 shows that whilst the global background correction gives 
a moderate probability for both the fit and the comparison 
with the composite, the local background correction method 
gives better agreement. This difference shows in the direct 
comparison of the magnitude distributions for the two back- 
ground correction methods (Table 4). This is clearly an ex- 
ample where global background correction gives rise to these 
deviations. We note that E394 is one of the poorest clusters 
in both our sample and that of Colless. Colless found that 
E124 was marginally discrepant when compared against the 
composite LF. We find no evidence for any difference. Lastly, 
E400, which we found in Section 5.2 to be potentially sig- 
nificantly different from the composite, Colless finds to be 
a reasonable match. However, since this cluster gives a dif- 
ference for both background subtraction methods, we find 
our result to be reliable. We also note that the smaller er- 
rors derived for by us for both E400 and El 24 allow us to 
make more definitive statements than Colless can for these 
clusters. 

Finally we note that unlike Colless, we find that all the 
subsets of the data we took are well fit by a one parameter 
Schechter LF, though they are not necessarily in agreement 
with the fit to the composite. 

6.2 Implications for cluster formation models 

There have been few detailed N-body simulations of the for- 
mation of clusters that have suitable dynamic range to map 
out the whole of the expected LF as a function of the cos- 
mological model assumed. For example, van Kampen (1995) 
makes simple predictions for the LF based on extensive sim- 
ulations, but is led to the conclusion that his modelling can- 
not correctly account for all of the merging that takes place 
within the cluster. His fit to the bright end of the luminosity 
function (which we might expect to be easier to model) is 
similar to the one observed. 

There have been studies that lead to simpler predic- 



tions as to the expected form of the LF according to the 
dominant processes that occur during cluster formation. The 
mode of cannabalism for the growth of cD galaxies as sug- 
gested by Hausman & Ostriker (1978) predicts that there 
should be differences between the luminosity function when 
broken down by Bautz-Morgan class. We see no real ev- 
idence for such a difference when we consider the range 
— 18 > M > —21. There is an indication for such a trend 
in the galaxies brighter than M = —21, in the sense that 
the Bautz-Morgan class I sample has fewer bright galaxies. 
However, better photometry of the bright galaxies would be 
required to confirm this result. Merritt (1983,1984) argued 
that most of the properties of clusters are essentially 'frozen 
in' when they form and that merging and tidal stripping 
processes are ineffective in rich cluster cores. This leads to 
the prediction that the LF should appear to be much the 
same from cluster to cluster. However, HST imaging data of 
moderate redshift clusters (e.g. Dressier et al. 1994) tends to 
rule against such a model, since it provides convincing evi- 
dence that mergers are important at least in the early stages 
of cluster evolution. One possibility is that the bulk of the 
change in the LF between clusters occurs only amongst the 
very faintest or brightest galaxies, the regime to which our 
data is insensitive. For example, it is known that some clus- 
ters show an upturn in the LF at faint magnitudes (e.g. 
Driver et al. 1994, Bernstein et al. 1995), and there has al- 
ways been marginal evidence for greater deviations in the LF 
when the brightest cluster galaxies are included (e.g. Lug- 
ger 1986 and the comments in Section 5.1), though we again 
caution against over-interpretation of results on the bright- 
est cluster members. However, our data clearly shows that 
in the region of the 'break' in the LF, the cluster-to-cluster 
variance is extremely small. 



7 CONCLUSIONS 

We have considered the cluster galaxy luminosity function 
as derived from data in the Edinburgh-Durham Southern 
Galaxy Catalogue and the Edinburgh-Milano Redshift Sur- 
vey. Our dataset is the largest multi-redshift cluster sample 
to be considered to date for this purpose. From this we find: 

(i) Our data shows that there is strong evidence that most 
clusters have very similar LFs when both the brightest and 
faintest galaxies are discarded. The average of the composite 
LFs derived using the 3 background correction methods in 
Section 5.1 give M* = -20.16 ± 0.02 and a = -1.22 ± 0.04, 
when the LF is determined in the range —18 > M > —21. 

(ii) We find a very good agreement between our results 
and those derived by Colless (1989) from similar input data, 
and with the mean of the data presented by Dressier (1978) 
and Lugger (1986). 

(iii) We have tested the stability of the derived Schechter 
luminosity function parameters to the type of background 
correction made. We find excellent agreement between a 
method based on a global background correction and that 
based on a strictly local correction. However, the largest 
systematic error remaining in our derivation of the best fit 
parameters for the composite LFs is the background correc- 
tion, with an uncertainty of up to 0.05 in M* and 0.05 in a 
possible. 
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(iv) We also tested the universality of the luminosity func- 
tion, by comparing individual cluster LFs with the compos- 
ite function derived in section 5.1, and by testing each clus- 
ter LF separately against a Schechter function (section 5.2). 
From this we find that at most ~10% of the clusters may 
have LFs that are significantly different from the composite. 

(v) We also broke our sample into sub-samples defined by 
richness, Bautz-Morgan class and velocity dispersion (sec- 
tion 5.3). From these tests we found weak evidence that 
the high and low velocity dispersion samples have different 
composite LFs, and that the higher velocity dispersion clus- 
ters may have a different LF to the global composite. This 
agrees with the marginal detection of a difference between 
similar samples by Colless (1989). We also found evidence 
for differences between the Bautz-Morgan class samples, but 
only for the brightest galaxies (M < —21) that we otherwise 
excluded for reasons of photometric reliability (Section 3.4). 
The trend evident is in line with the predictions of the canna- 
balism model of Hausman & Ostriker (1978), but we caution 
against over-interpretation of this result without more reli- 
able photometry for the bright galaxies. There is no evidence 
for any convincing difference between the other subsets of 
the data. Our results, and our simulations, do however show 
that it is important to consider more than one method of 
testing the difference between subsets of the data. On the 
basis of a one parameter fit alone we would have concluded 
that there were significant differences between the low and 
high velocity dispersion samples and the composite, and be- 
tween the low and high richness clusters. The other tests 
show no evidence for any difference however, with the ex- 
ception of the high velocity dispersion clusters where the 
two parameter fits, and Monte-Carlo simulations show weak 
evidence for a difference from the composite. 

It is likely that to make any further progress on the 
nature of the cluster galaxy luminosity function we will need 
to overcome the limitations of data such as ours. First, there 
is clearly a need to derive cluster membership free of any 
correction for backgrounds, by obtaining sufficient redshifts 
that we can actually map the galaxy distribution in three 
dimensions around the cluster. Secondly, better photometry 
is required to tie down the bright and faint ends of the LF. 
Given these it should be possible to not only discriminate 
between differences in the data broken down by richness, 
velocity dispersion or whatever else is desired, but also to 
test for the differences that our data suggest may exist in 
the LF when comparing the very brightest and very faintest 
galaxies. 
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Table 1. The full cluster sample. The EDCC number refers to the identification in Lumsdcn ct al. (1992). The field is the original 
Schmidt J survey field number. Ntot is the actual number of redshifts obtained for that cluster, and N c ; us is the number of those 
redshifts adjudged to lie in the cluster. BM class is the Bautz-Morgan class, where 1 represents class I, 3 represents class III, 2 represents 
any other defined class and implies that the cluster concerned is not in the Abell et al. (1989) catalogue and therefore has not been 
classified. The velocity dispersion is taken directly from Collins et al. (1995). 
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Table 2. The measured distribution of the fitted values of M* as a function of richness from Monte-Carlo simulations. The percentage 
values listed refer to the values of M that exclude that fraction of the actual measured fits. Since the distribution is skew, two limits 
are given for each percentage. A * implies that no acceptable fits could be found (M* < —26) for that fraction of the data. The last 
row gives the results for a simulation in which data was fitted over the range —23 to —18, showing the effect that the truncation at the 
bright end of the LF has on the fitting process. 
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Table 3. Fits to the various composite LFs described in the text. The columns are: the description of the composite; the derived values 
of M*and a (the error corresponds to the values of M* and a where \ 2 = X„jn + -*-)' tnc f° rma l likelihood that the fit is good; the 
number of clusters used in constructing the composite; the derived value of n*; and the measured richness of the composite. 
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Table 4. Best fit values of M* for a one parameter Schechter function fit to the individual cluster LFs. The error quoted corresponds to 
the mean of the deviations from M* where x 2 = xi,„ + 1- The probability quoted is the level at which this fit can be accepted. Both 
global and local background subtraction methods are given. For both of these the probability from a 2-sample x 2 tes t that the actual 
data differs from the composite LF is also given as p(CLF). Similarly, The probability that the local and global background subtracted 
data are the same is given as p(same). Lastly, the faint magnitude limit used in these tests is given. The bright limit is always —21. 
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Table 5. Fits to the subsets of the full cluster data as described in the text. The columns are: the description of the subset; the derived 
values of M*and a (again, the error corresponds to the values of M* and a where x 2 = X^in + •*■)> * ne f° rma l likelihood that the fit is 
good; the number of clusters that have been combined to form the composite for the subset; the derived value of n* ; the measured value 
of the cluster richness; lastly, the probability that the subset and the composite derived in section 4.1 have the same distribution using 
a 2-sample \ 2 test. The first row for each entry gives the best two parameter fit, the second gives a one parameter fit with a = —1.25. 
Those entries that appear only in the first row have the same values in the second. 
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Figure Captions 

Figure 1: The relationship between the catalogued EDSGC bj magnitude and the derived bj magnitude from the CCD 
sequences. The solid line represents the best fitting quadratic given in the text. The dashed line is given as a guide to the eye 
to show the deviation between the CCD magnitude scale and the EDSGC scale. 

Figure 2: The number counts from the EDSGC derived from the whole catalogue (the solid line), and from the average of all 
the annuli used for background correction of the LF (•). The difference between the two never exceeds 10%. The error bars 
shown represent the variance of the counts derived from the annuli. 

Figure 3: Derived values of M* from one parameter fits to model Schechter luminosity functions for clusters with (a) Ri = 10 
and (b) Ri = 50, when the fit is over the range -18 > M > -21. The model LF has M* = -20.2 and a = -1.25. When this 
range is extended to — 18 > M > —23, the distribution of M* for Ri — 25 is given in (c). This is clearly similar to the higher 
richness distribution when fitted over the narrower range in absolute magnitude. 

Figure 4: Contours of the derived distribution in M* and a for 10000 simulations of a two parameter fit to the standard 
Schechter function. The initial model LF has M* = -20.2 and a = -1.25. The contours contain 10%, 68.5%, 90%, 95% 
and 99% of the data. The initial richness given to the model LF is R c = 250, close to the measured value for most of our 
subsamples in section 5.3. This is therefore the correct error map for comparison with the observed values of M* and a from 
those composite sub-samples. 

Figure 5: Measured luminosity functions for (a) the local background correction method, (b) global background correction 
and (c) global correction but using the average of the annuli from (a). The data are fitted to a Schechter function in the 
range —21 < M < —18. The best fitting Schechter function is shown as a solid line. Only the data derived assuming a 50% 
background variance is shown. 

Figure 6: Measured luminosity functions for all those cluster with measured richness values, Ri, of 15 or more using the local 
background correction method. The best one parameter fit to the Schechter function is shown by a solid line. The fit is over 
the range —21 < M < M Um . 

Figure 7: Luminosity functions for (a) Bautz-Morgan class 1 clusters with Ri > 20, (b) Bautz-Morgan class 3 clusters with 
Ri > 20, (c) Bautz-Morgan class 1 clusters with Ri > 15, (d) Bautz-Morgan class 3 clusters with Ri > 15, (e) Ri > 20 clusters, 
(f) 10 < Ri < 20 clusters, (g) clusters with velocity dispersion < 700kms~ 1 and Ri > 20, (h) clusters with velocity dispersion 

> 700kms~ 1 and Ri > 20, (i) clusters with velocity dispersion < 700kms~ 1 and Ri > 15, (j) clusters with velocity dispersion 

> 700kms _1 and Ri > 15. The best fitting two parameter fits to a Schechter function over the range —18 > M > —21 are 
shown as solid lines. The best fitting one parameter fits (with a — —1.25) are shown as dashed lines. As can be seen, there is 
often a marked difference between the best fitting one and two parameter fits. 
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